Method for active noise cancellation using independent component analysis

ABSTRACT

The present invention considers a method for active cancellation using independent component analysis. More particularly, the present invention relates to a method which is operable the independent component analysis technique to an adaptive algorithm that can consider secondary or more higher statistical characteristics. 
     The conventional active noise cancellation systems mainly use the LMS(Least Mean Square) which considers secondary statistics among input signals. 
     Being different from the conventional active noise cancellation systems, the present invention provides a method for active noise cancellation using independent component analysis, which makes output signals independent of each other by considering secondary or more higher statistical characteristics. 
     Therefore, according to the present invention, the improved performances of the noise cancellation systems can be provided compared with the conventional active noise cancellation system which uses the LMS adaptive algorithm.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention considers a method for active noise cancellation using independent component analysis. More particularly, the present invention relates to a method which is operable the independent component analysis technique to an adaptive algorithm that can consider secondary or more higher statistical characteristics.

2. Description of the Related Art

FIG. 1 shows a structure of a general active noise cancellation system. In FIG. 1, a signal source(10) s is transmitted to a sensor through a channel, and a noise source(20) n₀ is input in the sensor so that the combined signal and noise s+n₀ form primary input(30) in the noise cancellation system.

The secondary sensor receives noise n₁ through another channel, and the sensor forms reference input(40) in the noise cancellation system. The noise n₁ is filtered to produce an output z, which is as close as possible of n₀, by passing through an adaptive filter(50), and the primary input s+n₀ deducts output z through an adder(60) and forms system output(70) u=s+n₀−z in the noise cancellation system.

The purpose of the conventional active noise cancellation is to get output u=s+n₀−z which is as close as possible of signal s in the point of least squares. To reach the purpose, the filter is adapted using a least mean square(LMS) adaptive algorithm to minimize the entire output of the noise cancellation system. In other words, the output in the active noise cancellation system is operated as an error signal during adaptation.

The coefficient adaptation of the filter follows the Widrow-Hoff LMS algorithm and can be expressed as following expression. Δw(k)αu(t)n ₁(t−k)  [Expression 1]

Where, w(k) is k^(th) order coefficient, and t is sample index.

SUMMARY OF THE INVENTION

The purpose of the present invention is to provide a method for active noise cancellation using independent component analysis, which allows to get the improved performances of the noise cancellation systems compared with the conventional LMS adaptive algorithm.

To reach the said purpose of the present invention, being different from the conventional noise cancellation systems, the present invention provides a method for active noise cancellation using independent component analysis, which makes output signals independent of each other by considering secondary or more higher statistical characteristics.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 shows a structure of a general active noise cancellation system.

FIG. 2 shows a structure of a feedback filter according to the present invention.

FIG. 3 shows a structure of a feedback filter for active noise cancellation according to the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

Hereafter, an embodiment according to the present invention is described in detail by referring to accompanying drawings.

FIG. 2 shows a structure of a feedback filter according to the present invention, and FIG. 3 shows a structure of a feedback filter for active noise cancellation according to the present invention.

Independent component analysis can recover unobserved independent source signals from the combined signals of sound sources which are mixed through an unknown channel. Let us consider a set of unknown sources, s(t)=[s₁(t), s₂(t), . . . , s_(N)(t)]^(T), such that the components s_(i)(t) are zero-mean and mutually independent. From N sensors, it is obtained a set of signals, x(t)=[x₁(t), x₂(t), . . . , x_(N)(t)]^(T), which are mixed through a channel from the sources, s(t)=[s₁(t), s₂(t), . . . , s_(N)(t)]^(T). If a channel can be modeled as instantaneous mixing, the sensor signals, x(t)=[x₁(t), x₂ (t), . . . , x_(N)(t)]^(T), can be expressed as following expression. x(t)=A·s(t)  [Expression 2]

Where, A is an unknown invertable matrix called as a mixing matrix.

Therefore, it is the subject to retrieve the sound sources by finding the inverse matrix of the mixing matrix A using only sensor signals x(t).

However, recovering a permuted and rescaled version of original sound sources doesn't matter from the standpoint of source separation because it doesn't affect their waveforms.

Therefore, by estimating the unmixing matrix W, the retrieved signals u(t), which are the original ones up to permutation and scaling, is acquired from the following expression. u(t)=W·x(t)  [Expression 3]

Herein, to estimate the unmixing matrix W, it is assumed that each of sound sources is independent. It means that a signal from one sound source doesn't affect a signal from another source, and it can be considered as possible suppose. Also, the statistical independence includes the statistical characteristics of all orders.

Because the statistical independence doesn't affect permutation and scaling, the unmixing matrix W of expression 3 can be obtained. The unmixing matrix W is learned by using the following expression which makes the statistical independence among estimated source signals maximize. $\begin{matrix} {{{\Delta\; W} \propto {\left\lbrack W^{T} \right\rbrack^{- 1} - {{\varphi(u)}x^{T}}}},{{\varphi\left( {u_{i}(t)} \right)} = {- \frac{\frac{\partial{P\left( {u_{i}(t)} \right)}}{\partial{u_{i}(t)}}}{P\left( {u_{i}(t)} \right)}}}} & \left\lbrack {{Expression}\mspace{14mu} 4} \right\rbrack \end{matrix}$

Where, P(u_(i)(t)) approximates the probability density function of estimated source signal u_(i)(t).

In real world situations, instantaneous mixing is hardly encountered, and the mixing of sources involves convolution and time-delays as following expression. $\begin{matrix} {{x_{1}(t)} = {\sum\limits_{j = 1}^{N}\;{\sum\limits_{k = 0}^{K - 1}\;{{a_{ij}(k)}{s_{j}\left( {t - k} \right)}}}}} & \left\lbrack {{Expression}\mspace{14mu} 5} \right\rbrack \end{matrix}$

Where, x_(i)(t) is a measured sensor signal, s_(j)(t) is a sound source signal, and a_(ij)(k) is a coefficient of a mixing filter of length K.

To separate the source signals, the feedback filter structure of FIG. 2 can be used. Where, the recovered signal u_(i)(t) can be expressed as following. $\begin{matrix} {{u_{i}(t)} = {{\sum\limits_{k = 0}^{K}\;{{w_{ij}(k)}{x_{1}\left( {t - k} \right)}}} + {\sum\limits_{j = {{1j} \neq i}}^{N}\;{\sum\limits_{k = 1}^{K}\;{{w_{ij}(k)}{u_{j}\left( {t - k} \right)}}}}}} & \left\lbrack {{Expression}\mspace{14mu} 6} \right\rbrack \end{matrix}$

Where, w_(ij)(k) shows k^(th) order coefficient of the filter to estimate an original sound sources.

In this structure, there are three different cases of the filter coefficients. In other words, these are w_(ii)(0) which is a zero delay coefficient in a direct filter, w_(ii)(k), k≠0 which is a delay coefficient in a direct filter, and w_(ij)(k), i≠j which is a coefficient in a feedback cross filter. Learning rules for all these cases are as following. $\begin{matrix} {\mspace{79mu}{{{\Delta\;{w_{ii}(0)}} \propto {{1/{w_{ii}(0)}} - {{\varphi\left( {u_{i}(t)} \right)}{x_{i}(t)}}}},\mspace{79mu}{{\Delta\;{w_{ii}(k)}} \propto {{- {\varphi\left( {u_{i}(t)} \right)}}{x_{i}\left( {t - k} \right)}}},\mspace{79mu}{{\Delta\;{w_{ij}(k)}} \propto {{- {\varphi\left( {u_{i}(t)} \right)}}{u_{j}\left( {t - k} \right)}}},{{\varphi\left( {u_{i}(t)} \right)} = {- \frac{\frac{\partial{P\left( {u_{i}(t)} \right)}}{\partial{u_{i}(t)}}}{P\left( {u_{i}(t)} \right)}}}}} & \left\lbrack {{Expression}\mspace{14mu} 7} \right\rbrack \end{matrix}$

Zero delay coefficient, w_(ii)(0), scales the data to maximize the information transmitted through the nonlinear function, delay coefficient, w_(ii)(k), k≠0, whitens each output from the corresponding input signal temporally, and coefficient in a feedback cross filter, w_(ij)(k), i≠j, decorrelates each output φ(u_(j)(t)) from all other recovered signal u_(j)(t).

FIG. 2 shows the feedback filter structure with two inputs and outputs. However, by increasing the number of inputs and outputs in the feedback filter structure, a feedback filter structure can be constructed with an arbitrary number of inputs and outputs. And, in this structure, the same expressions as expression 5, 6, and 7 that show the mixed signals, recovered signals, and learning rules of filter coefficients can be used.

For the convolved mixtures, the recovered signals are the whitened origin source signals temporally. Therefore, whitening problem can be solved by remaining zero delay coefficient, w_(ii)(0), and fixing delay coefficient, w_(ii)(k), k≠0.

In case that pure noise can be observed, noise component from mixed signal can be canceled by using independent component analysis.

In FIG. 1, which shows basic active noise cancellation, it can be assumed that noise n₁ and signal s are independent to each other, but noise n₁ is related to noise n₀ in unknown way. Because signal and noise do not affect each other, this assumption is reasonable.

Independent component analysis can be used to cancel noise dependent component of primary input by the secondary input. The primary input is the mixed signal of noise and signal, and the secondary input is noise component without signal component, therefore, the feedback filter structure for independent component analysis can be modified.

In other words, when the primary and secondary input are x₁ and x₂ respectively like in FIG. 2, the secondary input x₂ is the noise component without signal component, so it is not required feedback filter W₂₁(k) of expression 6.

The feedback filter structure, which has two inputs and outputs that do not use feedback filter w₂₁(k), is shown as the structure of FIG. 3 and can be expressed as following expression. $\begin{matrix} {{{u_{1}(t)} = {{\sum\limits_{k = 0}^{K}\;{{w_{11}(k)}{x_{1}\left( {t - k} \right)}}} + {\sum\limits_{k = 1}^{K}\;{{w_{12}(k)}{u_{2}\left( {t - k} \right)}}}}},{{u_{2}(t)} = {\sum\limits_{k = 0}^{K}\;{{w_{22}(k)}{x_{2}\left( {t - k} \right)}}}}} & \left\lbrack {{Expression}\mspace{14mu} 8} \right\rbrack \end{matrix}$

By comparing FIG. 1 which shows basic active noise cancellation system and FIG. 3, the mixture s+n₀ of signal and noise, which forms the primary input, is corresponding to x₁, whereas n₁, which forms the secondary input, is corresponding to x₂. Also, output u is corresponding to u₁.

Herein, in case of remaining zero delay coefficient, w_(ii)(0), and fixing delay coefficient, w_(ii)(k), k≠0to 0, because feedback filter w₂₁(k) is not used, output signal u₁ of expression 8 corresponding to output u of the conventional active noise cancellation system can be expressed using two inputs x₁ and x₂ as following. $\begin{matrix} {{u_{1}(t)} = \;{{{w_{11}(0)}{x_{1}(t)}} + {{w_{22}(0)}{\sum\limits_{k = 1}^{K}\;{{w_{12}(k)}{x_{2}\left( {t - k} \right)}}}}}} & \left\lbrack {{Expression}\mspace{14mu} 9} \right\rbrack \end{matrix}$

The second term of the expression 9 corresponds to output z of the adaptive filter in the conventional active noise cancellation system. Therefore, the feedback filter structure of independent component analysis for active noise cancellation can be considered as the same structure of the conventional active noise cancellation system.

The learning rule of each filter coefficient for an active noise cancellation system can be expressed by independent component analysis as following. $\begin{matrix} {{{\Delta\;{w_{u}(O)}} \propto {{1/{w_{u}(O)}} - {{\varphi\left( {u_{i}(t)} \right)}{x_{i}(f)}}}},} & \left\lbrack {{Expression}\mspace{14mu} 10} \right\rbrack \\ {{{\Delta\;{w_{u}(k)}} \propto {{- {\varphi\left( {u_{i}(r)} \right)}}{x_{i}\left( {t - k} \right)}}},} & \; \\ {{{\Delta\;{w_{ij}(k)}} \propto {{- {\varphi\left( {u_{i}(t)} \right)}}{u_{j}\left( {t - k} \right)}}},{{\varphi\left( {u_{i}(t)} \right)} = \frac{\frac{\partial{P\left( {u_{i}(t)} \right)}}{\partial{u_{i}(t)}}}{P\left( {u_{i}(t)} \right)}}} & \; \end{matrix}$

Output of the active noise cancellation system obtained by learning filter coefficients is signal component which is independent of noise component.

In FIG. 3, the feedback filter structure for active noise cancellation with two inputs and outputs can be extended to an structure with arbitrary number of inputs and outputs like in FIG. 2, and the recovered signals and learning rules of coefficients in the said feedback filter structure for active noise cancellation which has arbitrary number of inputs and outputs can be obtained by extending expression 8 and 10 in the same way.

As mentioned above, the method for active noise cancellation using independent component analysis according to the present invention provides the improved noise cancellation performances compared with those of the active noise cancellation system which uses the conventional LMS adaptive algorithm. 

1. The method for active noise cancellation using independent component analysis which is characterized by controlling active noise by learning each adaptive filter coefficient according to the following expression in active noise cancellation system, wherein the mixture x_(i) of signal and noise that forms the primary input and noise x₂ that forms the secondary input is related by the following expression: $\begin{matrix} {\mspace{79mu}{{{\Delta\;{w_{ii}(0)}} \propto {{1/{w_{ii}(0)}} - {{\varphi\left( {u_{i}(t)} \right)}{x_{i}(t)}}}},\mspace{79mu}{{\Delta\;{w_{ii}(k)}} \propto {{- {\varphi\left( {u_{i}(t)} \right)}}{x_{i}\left( {t - k} \right)}}},\mspace{79mu}{{\Delta\;{w_{ij}(k)}} \propto {{- {\varphi\left( {u_{i}(t)} \right)}}{u_{j}\left( {t - k} \right)}}},{{\varphi\left( {u_{i}(t)} \right)} = {- \frac{\frac{\partial{P\left( {u_{i}(t)} \right)}}{\partial{u_{i}(t)}}}{P\left( {u_{i}(t)} \right)}}}}} & \left\lbrack {{Expression}\mspace{14mu} 7} \right\rbrack \end{matrix}$ wherein, the said w_(ii)(0) is a zero delay coefficient in a direct filter, w_(ii)(k), k≠0 is a delay coefficient in a direct filter, w_(ij)(k), i≠j is a coefficient in a feedback cross filter, Δ before each coefficient is a change amount of the corresponding coefficient, t is sample index, and P(u_(i)(t)) approximates the probability density function of estimated source signal u_(i)(t).
 2. The method for active noise cancellation using independent component analysis according to claim 1, which is characterized by obtaining the said u_(i)(t) by following expression: ${{u_{1}(t)} = {{\sum\limits_{k = 0}^{K}\;{{w_{11}(k)}{x_{1}\left( {t - k} \right)}}} + {\sum\limits_{k = 1}^{K}\;{{w_{12}(k)}{u_{2}\left( {t - k} \right)}}}}},{{u_{2}(t)} = {\sum\limits_{k = 0}^{K}\;{{w_{22}(k)}{{x_{2}\left( {t - k} \right)}.}}}}$
 3. The method for active noise cancellation using independent component analysis according to claim 1, wherein signal cancellation range corresponding to active noise is extended for the system which acquires many noise signals or mixtures of signal and noise by increasing the number of inputs or outputs of the said active noise cancellation system. 